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EFFECTS OF WEAK LENSING ON THE TOPOLOGY OF CMB MAPS 
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ABSTRACT 

We investigate the non-Gaussian signatures in Cosmic Microwave Background (CMB) maps induced 
by the intervening large-scale structure through weak lensing. In order to measure the deviation from 
the Gaussian behavior of the intrinsic temperature anisotropics, we use a family of three morphological 
descriptors, the so-called Minkowski functionals. We show analytically how these quantities depend 
on the temperature threshold, and compare the results to numerical experiments including the instru- 
mental effects of Planck. Minkowski functionals can directly measure the statistical properties of the 
displacement field and hence provide useful constraints on large-scale structure formation in the past. 

Subject headings: Methods: statistical — cosmic microwave background — dark matter — 
gravitational lensing — large-scale structure of Universe 



1. INTRODUCTION 

The temperature anisotropy in the CMB is a powerful 
probe of the content and nature of our Universe. Most 
inflationary scenarios (Sato 1981; Guth 1981) predict that 
the temperature fluctuation field obeys Gaussian statistics 
(Guth & Pi 1985), so all its statistical properties can be ac- 
curately predicted based on the Gaussian theory (Bardeen 
et al. 1986; Bond & Efstathiou 1987). However, gravita- 
tional lensing by the matter inhomogeneities between the 
last scattering surface and us imprints non-Gaussian sig- 
natures on the CMB. These signatures directly probe the 
mass distribution up to very high redshift, and their mea- 
surement could greatly help to constrain cosmological pa- 
rameters. 

Various statistical methods have been used to investi- 
gate lensing of CMB maps. The effect on the power spec- 
trum Ce itself is rather small (Seljak 1996). The prob- 
ability density function (PDF) of peak ellipticities of the 
lensed temperature field is in principle sensitive to the lens- 
ing signatures on scales below 10', but the finite beam size 
of detectors tends to circularize the deformed ellipticities 
again (Bernardeau 1998). Using the correlation between 
ellipticities of the lensed CMB map and distant galaxies 
proves more robust against the beam smearing effect (van 
Waerbeke et al. 2000). Takada et al. (2000) measure lens- 
ing signatures on scales around 75' with the two-point cor- 
relation function of hotspots. 

In this Letter, we suggest to detect weak lensing sig- 
natures in CMB maps with the Minkowski functionals 
(Minkowski 1903). To illustrate that the method is useful 
at all, we perform quantitative analyses on numerical ex- 
periments (Takada & Futamase 2000). Section 2 describes 
the production of the CMB maps used in our work. In 
Section 3, we summarize some of the properties of the 
Minkowski functionals and apply them to the maps. Sec- 
tion 4 provides an outlook. The Appendix briefly derives 
the dependence of the Minkowski functionals of a lensed 



CMB map on the temperature threshold. 

2. NUMERICAL EXPERIMENTS 

We consider two currently fashionable cosmologies, 
namely Cold Dark Matter variants with f2 m =l and h=0.5 
(SCDM) and Sl m =0.3, n x =0.7, and h=0.7 (ACDM). Q m 
and fl\ are the energy densities of matter and vacuum 
relative to the critical density, respectively, and h is the 
Hubble constant in units of 100km s _1 Mpc _1 . The baryon 
density is chosen as Sl\,h 2 = 0.019 from recent Big Bang 
Nucleosynthesis results (Buries & Tytler 1998). We use a 
Harrison-Zel'dovich initial power spectrum and the trans- 
fer function of Bardeen et al. (1986) with the shape param- 
eter form (Sugiyama 1995). After that, we can only vary 
the normalization, or equivalently as, the mass fluctuation 
in a sphere of radius 8/i^ 1 Mpc. 

For each model, we calculate the Cg with CMBFAST 
(Seljak & Zaldarriaga 1996) and compute 30 realizations 
of the CMB on an 80° x 80° square divided into 2' pixels. 

The lensing effect can be treated as a mapping of the 
primary CMB anisotropies. When we observe the CMB 
temperature tp in a certain direction 6 on the sky, we ac- 
tually see the intrinsic field (p at a different position 6 + £ 
on the last scattering surface: 

iK0) = <p(6+t(6))' (!) 
To obtain the lensing displacement field £, we first gener- 
ate the convergence field k. Its power spectrum P K (£) is 
the matter power spectrum Pg{k) projected along the line 
of sight: 



P K {1) = \Hfo r 
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X and r(x) denote the comoving radial distance and the 
corresponding comoving angular diameter distance, and 
Xrcc is the radial distance to the last scattering surface. 
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From k we compute the displacement field £ in Fourier 
space through £(k) = 2ikfc~ 2 ft(k) and transform to real 
space. According to Equation (1), the displacement field 
only distorts the locations of the values in the primordial 
map. So we interpolate the original temperature values 
with a cloud-in-cell procedure to obtain the lensed temper- 
ature map on a regular grid. Finally, we mimic an obser- 
vation through Planck (Mandolesi et al. 1995) by adding 
a Gaussian beam smearing and white noise. The beam 
FWHM and the noise level per FHWM pixel are chosen 
as 5.5' and 4.3 x 1CT 6 , respectively (Efstathiou & Bond 
1999). 

As an example, Figure 1 shows a map with and without 
lcnsing. 

3. INTEGRAL GEOMETRY FOR CMB MAPS 

Minkowski functional were formally introduced into 
cosmology by Mecke et al. (1994). Although the whole 
family of these geometric descriptors was already men- 
tioned in (Coles 1988), Schmalzing & Gorski (1998) first 
applied them to CMB maps. Bernardeau (1997) sug- 
gested Minkowski functional for detecting weak lcnsing 
signatures in the CMB. For a continuous field, we use the 
Minkowski functionals to describe the geometry of its ex- 
cursion sets, that is the area above a given temperature 
threshold. In two dimensions, there are three Minkowski 
functionals which correspond to well-known geometrical 
quantities, namely the area i>o, the circumference vi, and 
the Euler characteristic v 2 - 

The unlensed CMB maps obey Gaussian statistics. 
Therefore, their statistical properties can be completely 
characterized by the power spectrum Cg. Following 
Tomita (1990), we parametrize the Minkowski functionals 
of the isotemperature contour at the normalized threshold 
v with r, the variance of the field's first derivative 3 
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The Appendix shows that the shapes of the Minkowski 
functional curves of the lensed field remain the same as in 
the Gaussian case. So the non-Gaussianity manifests it- 
self only in the normalization of the curves. Let us denote 
the constants of proportionality from Equations (A4) and 
(A6) with «i and a 2 , respectively. Equation (3) yields 



Gauss 



, Gauss 



(4) 



Therefore, for the Gaussian random field the quantity 

vw a 2 



(5) 



equals unity. As we shall see, lcnsing leads to deviations 
from unity. 

Even with high resolution and large sky coverage, these 
changes are easily overcome by the cosmic error. There- 
fore, accurate evaluation of the Minkowski functionals is 
crucial. We use three independent methods and take 
care that they produce compatible results. One of them 
uses contour integration to evaluate the circumference and 
the Morse theorem (Morse & Cairns 1969) to determine 
the Euler characteristic (see e.g. Novikov et al. 1999). 
Two more elaborate methods (Schmalzing & Buchert 
1997) evaluate the Minkowski functionals via Crofton's 
formula (Crofton 1868), and by averaging over invariants 
formed from derivatives, respectively. Figure 2 shows the 
Minkowski functionals of one of the models. 

For each realization, we determine the non-Gaussianity 
parameter a by fitting the measured Minkowski function- 
als to their expected Gaussian shapes. Since each of our re- 
alizations covers only 15% of the sky, while for the Planck 
satellite a sky coverage of 60% and more is expected, we 
reduce our estimated variances by a factor of two. It is 
also worth mentioning that we assume Gaussianity of the 
convergence field on scales of 5' and above. This is not ex- 
actly true (Jain et al. 2000), so we expect an even stronger 
signal from a refined analysis with a convergence field cal- 
culated from A-body simulations. 

Table 1 summarizes our results. As expected, neither of 
the unlensed models deviates significantly from Gaussian- 
ity. For the maps that include the weak lensing effect, the 
average a is different from one. In all but one case, this 
difference is significant, and in two out of the four investi- 
gated models with lensing, the significance level is above 
95%. 

4. DISCUSSION AND OUTLOOK 

We have measured the weak lensing effect of large-scale 
structure on the observed temperature anisotropies of the 
CMB with Minkowski functionals. Numerical simulations 
have shown that the effect can be significant when ob- 
served with the experimental specifications of Planck. It 
remains to be seen whether the Minkowski functionals can 
directly measure any characteristics of large-scale struc- 
ture. Since they are sensitive to smoothing, we expect that 
varying the smoothing scale can reveal information on the 
convergence field on different scales. Most importantly, 
however, we have proven that this method can measure 
non-Gaussian signatures induced by weak lensing at all. 
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APPENDIX 
ANALYTICAL TREATMENT 

Consider a scalar random field -0 in two dimensions, e.g. the observed CMB temperature anisotropy. It can be related 
to a scalar Gaussian random field <p, the intrinsic temperature anisotropy, through the vector-valued random field the 
displacement, by Equation (1). 

It is well-known that the PDF of the temperature field, and hence the zeroth Minkowski functional vq, which is just the 
integrated PDF, does not change at all under lensing. We write the other two Minkowski functionals as spatial averages 
over invariants formed from the field's derivatives 4 (Schmalzing & Gorski 1998): 

t\ n /xr, n / ,2 , ,2 \ / \ 1 I , 2^,1^,2^,12 - ^1^,22 - ^ 2 2 V>,n \ , A1 s 
«iW = ^^-f)^ 2 i+^), v 2 {v) = — ( 5{$-v) + ^ 2 : ). (Al) 

In order to evaluate these two averages, we need to express the first- and second-order derivatives of -0 in terms of the 
fields if and £. Straightforward differentiation yields 5 

= </?,i + f,k£,k,i, 4>,ij = <fi,ij + ¥dk£,k,j + ¥>,jk£k,i + <P,kl£k,i£l,j + <P,k£,k,ij- (A2) 

For the circumference v\, Equation (Al) involves only the first derivatives of the observed field ip and therefore, by 
Equation (A2), only the first derivatives of the intrinsic field ip. Since these are independent of the value <p itself (Adler 
1981), and £ is of course independent of the field ip, the average splits neatly into two factors: 

«iM - \ W - *)) + - J 2 -^ ( V%+^> • ( A3 ) 

The remaining average is independent of v, so the curve has the Gaussian shape: 

vi (i/) oce"" 2 / 2 . (A4) 

Turning to the Euler characteristic we observe that Equation (Al) expressed in terms of the fields <p and ^ depends 
linearly on the second derivatives tp^j and £k,ij ■ Therefore, v 2 depends on the threshold v only through 

e — 2 /2 

(S((f - v)<p,ij) = -tv — j^Sij and (S(ip - v)£k,ij) = 0. (A5) 



2tt 



The remainder of the average in Equation (Al) again produces factors that do not depend on the threshold v. So the 
curve v<2.(y) also has the shape expected for the Gaussian case: 

v 2 {y) oc vc-" 2 ' 2 . (A6) 



4 Indices following a comma denote a spatial derivative. 
5 Summation over pairwise indices is implied throughout. 
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Table 1 

The "non-Gaussianity parameter" a from Equation (5). For a Gaussian random field, a should be unity, while 

DEVIATIONS FROM ONE ARE EXPECTED FOR A WEAKLY LENSED CMB SKY. a AND Aa ARE THE AVERAGE AND STANDARD ERROR 

OVER ALL REALIZATIONS OF EACH MODEL. 
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Fig. 2. — Minkowski functionals for the lensed and unlensed maps observed by the Planck satellite. Results for the unlensed maps are 
shown in grey, while black lines indicate the Minkowski functionals of the lensed maps. For both, the average and standard deviation are 
shown in solid lines, while the expectation values for a Gaussian random field with the same two-point characteristics are displayed in dashed 
lines. Note that we only show the region around the maxima of the Minkowski functional curves, where the differences are visibly significant. 



